arXiv:1502.05399v2 [astro-ph.GA] 3 May 2015 


Draft version May 5, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


A SPECTROSCOPIC REDSHIFT MEASUREMENT FOR A LUMINOUS LYMAN BREAK GALAXY 

AT Z = 7.730 USING KECK/MOSFIRE 

P. A. Oesch^’^, P. G. van Dokkum^, G. D. Illingworth^, R. J. Bouwens^, I. Momcheva^, B. Holden^, 
G. W. Roberts-BorsanP’®, R. Smit®, M. Frank'*, I. Labbe*, V. Gonzalez^, D. Magee^ 

Draft version May 5, 2015 

ABSTRACT 

We present a spectroscopic redshift measurement of a very bright Lyman break galaxy at 2: = 
7.7302 ± 0.0006 using Keck/MOSFIRE. The source was pre-selected photometrically in the EGS field 
as a robust 2 ; ^ 8 candidate with H = 25.0 mag based on optical non-detections and a very red 
Spitzer/IRAC [3.6] —[4.5] broad-band color driven by high equivalent width [O IIl]+H/3 line emission. 
The Lyct line is reliably detected at 6.la and shows an asymmetric profile as expected for a galaxy 
embedded in a relatively neutral inter-galactic medium near the Planck peak of cosmic reionization. 
The line has a rest-frame equivalent width of EWq = 21 ±4 A and is extended with Efwhm = 36 OI 7 Q 
km s“^. The source is perhaps the brightest and most massive 2: ^ 8 Lyman break galaxy in the 
full CANDELS and BoRG/HIPPIES surveys, having assembled already Mq of stars at only 

650 Myr after the Big Bang. The spectroscopic redshift measurement sets a new redshift record 
for galaxies. This enables reliable constraints on the stellar mass, star-formation rate, formation 
epoch, as well as combined [O IIl]+H/d line equivalent widths. The redshift confirms that the IRAC 
[4.5] photometry is very likely dominated by line emission with EWo{[0 IIl]+H;d)= 7201^50 A. This 
detection thus adds to the evidence that extreme rest-frame optical emission lines are a ubiquitous 
feature of early galaxies promising very efficient spectroscopic follow-up in the future with infrared 
spectroscopy using JWST and, later, ELTs. 

Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: evolution — dark ages, 
reionization, first stars 


1 . INTRODUGTION 

The spectroscopic confirmation and characterization of 
galaxy candidates within the cosmic reionization epoch 
has been a major challenge for observational extragalac- 
tic astronomy for the last few years. Recently, large 
samples of several hundred galaxy candidates have been 
identified at 2 : ^ 7 — 11 thanks to the exceptional near- 
infrared sensitivity of the WEC3/IR camera onboard 
the Hubble Space Telescope {HST; e.g., Bouwens et al. 
2011, 2015; Schenker et al. 2013; McLure et al. 2013; 
Oesch et al. 2012 , 2014; Einkelstein et al. 2014). How¬ 
ever, despite this unprecedented target sample, very lit¬ 
tle progress has been made in spectroscopically confirm¬ 
ing galaxies in the cosmic reionization epoch. Gurrently, 
only a handful of normal galaxies have reliably measured 
redshifts at ^ > 7 (see, e.g., Vanzella et al. 2011; Pen- 
tericci et al. 2011 ; Ono et al. 2012 ; Schenker et al. 2012 ; 
Shibuya et al. 2012 ; Einkelstein et al. 2013), with most 
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spectroscopic surveys being unsuccessful or only result¬ 
ing in uncertain candidate lines (e.g. Treu et al. 2013; 
Jiang et al. 2013a; Tilvi et al. 2014; Garuana et al. 2014; 
Eaisst et al. 2014; Vanzella et al. 2014; Schenker et al. 
2014). 

Gurrent studies at ^ > 6 thus rely on photometric 
samples, with selection criteria and photometric redshifts 
that are somewhat uncertain due to the difficulty of es¬ 
tablishing reliable priors of potential contaminant pop¬ 
ulations at lower redshift. Spectroscopic follow-up is 
therefore particularly important for the very rare galaxies 
at the bright end of the UV luminosity and mass func¬ 
tions where any contamination has a very large impact. 

The low success rate of spectroscopic follow-up sur¬ 
veys is likely caused by a decreased fraction of galaxies 
showing strong Lya emission due to an increased neutral 
fraction in the inter-galactic medium (IGM) at 2 : > 6 
(Stark et al. 2010, 2011; Treu et al. 2013; Schenker et al. 
2012 , 2014; Pentericci et al. 2014). While Lya is the pri¬ 
mary spectral feature for spectroscopic confirmation of 
high-redshift candidates, new surveys targeting the Ly¬ 
man continuum break are underway using the WEG3/IR 
grism on the HST, or alternatively, weak UV lines may 
also be detectable from the ground (Stark et al. 2014). 

Two recent successful Lya detections of Lyman break 
selected galaxies at ^ > 7.2 were published in Ono et al. 
( 2012 , z = 7.213) and Einkelstein et al. (2013, ^ = 7.508). 
Both these sources are relatively bright with H = 25.2, 
and 25.6 mag and both show a significant fiux excess 
in their IRAC photometry, which is consistent with ex¬ 
tremely strong [O HI] line emission at 4.1 — 4.3 jam. Such 
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Fig. 1.— Top - Images showing a 5"x5" region around our 
primary target galaxy EGS-zs8-l in the HST and Spitzer/IRAC 
filters. These are, from left to right, Veoe, -1^814, Tios, J 125 , 
and 3.6 pLvo and 4.5 pirn.. Bottom - The spectral energy distri¬ 
bution of EGS-zs8-l based on fits to the HST+Spitzer+K-band 
photometry. Downward pointing vectors represent 2 cr upper lim¬ 
its in non-detection bands. A significant flux excess in the IRAC 
4.5 pivo band is evident. Together with the strong spectral break, 
this constrains the photometric redshift to ZpPot = 7.7 ± 0.3, in 
excellent agreement with the spectroscopic measurement as shown 
later. The best fit low-redshift solution at 2 : ~ 1.8 is shown as a 
gray line for completeness. However, this SED has a likelihood of 
<10“^ and is ruled out by the photometry. 


strong lines are characteristic of early star-forming galax¬ 
ies, as evidenced by a clear increase in broad-band flux 
excess with redshift (e.g., Schaerer & de Barros 2009; 
Labbe et al. 2013; Stark et al. 2013; Gonzalez et al. 
2014; Smit et al. 2014). Such excesses can be used to 
select relatively clean samples of star-forming galaxies at 
2 ; ^ 6.6 — 6.9 as well as 2 ; ^ 8 (e.g. Smit et al. 2015, 
Roberts-Borsani et ah, 2015, in prep.). 

In this Letter we present a successful spectroscopic red¬ 
shift measurement at 2 : ^ 8 using Keck/MOSFIRE of one 
of the brightest Lyman Break galaxies (LBGs) at that 
epoch. This galaxy was pre-selected as a high-priority 
target because of a very red [3.6]-[4.5] color, likely caused 
by strong [O III] emission. Our target selection is sum¬ 
marized in Section 2, while Section 3 outlines the spec¬ 
troscopic observations, and our results are presented in 
Section 4. 

Throughout this paper, we adopt LIm = 0.3, LIa = 
0.7, iLo = 70 kms“^Mpc“^, i.e. h = 0.7, consistent 
with the measurements from Planck (Planck Gollabora- 
tion et al. 2015). Magnitudes are given in the AB system. 

2. TARGET SELEGTION 

We briefly summarize our selection of a robust ^ ^ 8 
LBG sample over the GANDELS fields using extreme 
IRAG photometry. Eor more details see Roberts-Borsani 
et al. (2015, in prep.). 

The selection builds on Smit et al. (2015), who iden¬ 
tify a sample of 2 ; ^ 6.8 galaxies based on strong 
[O III] AA4959, 5007 plus H/3 emission lines resulting in 
very blue [3.6] —[4.5] IRAG colors. As these lines shift 
into the IRAG 4.5 /im band, galaxies at 2 : ^ 7 to 2 : ^ 9 
exhibit red [3.6] —[4.5] IRAG colors (see also Stark et al. 
2013; Labbe et al. 2013; Bowler et al. 2014). 


We exploit the availability of deep Spitzer/IRAG pho¬ 
tometry over the HST GANDELS-Wide fields to system¬ 
atically search for bright galaxies with IRAG colors of 
[3.6] —[4.5]> 0.5 mag in addition to a hya break (i.e., a 
non-detection at < 1 /im), characteristic for ^ > 7 galax¬ 
ies. This resulted in two candidates with H < 25.1 mag 
in the EGS field (see Roberts-Borsani et ah, in prep.). 
Eortuitously, these two sources are < 6 ' from each other 
and can be targeted in a single MOSEIRE mask. 

Stamps and SED fits for one of these sources (EGS- 
zs 8 -l) are shown in Eigure 1. The E606W, E814W, 
E125W and E160W images come from the GANDELS 
survey (Grogin et ai. 2011 ), while the IRAG images are 
from the SEDS survey (Ashby et al. 2013). Also shown 
are WEG3/IR E105W observations that are fortuitously 
available over this source as a result of a separate follow¬ 
up program (GO: 13792, PI: Bouwens). Despite its mod¬ 
est depth, this Yios image still provides a highly improved 
photometric redshift measurement by constraining the 
spectral break at 1 /im. 

As seen from the figure, the source EGS-zs 8 -l is only 
detected at > 1 /im in the WEG3/IR imaging as well as in 
both IRAG 3.6 and 4.5 /im bands. The [3.6] —[4.5] color of 
this source is measured to be 0.53T0.09, i.e., at the edge 
of our IRAG color selection window ([3.6] —[4.5] > 0.5). 

3. OBSERVATIONS 

3.1. MOSFIRE Spectroscopy 

We use the Multi-Object Spectrometer for Infra-Red 
Exploration (MOSEIRE; McLean et al. 2012) on the 
Keck 1 telescope for Y-band spectroscopy of our pri¬ 
mary 2: ^ 8 targets in the search for their Lya emission 
lines. MOSEIRE offers efficient muitipiex observations 
over a fieid of view of ^ 6 ' x 3' at a spectrai resoiution 
of R - 3000 {R = 3500 and R = 2850 for a 0'.'7 or 0'.'9 
siit, respectiveiy). 

Data over the EGS fieid were taken during three nights, 
2014 April 18, 23, and 25. While the first night was es- 
sentiaiiy iost due to bad seeing and ciouds, the remain¬ 
ing nights had better conditions with a median seeing of 
1 " and oniy few cirrus ciouds during the iast night. A 
dome shutter break probiem aiso ied to some vignetting 
during the iast 30 minutes of the April 25 night before 
we stopped observations early. In total, we obtained 2.0 
hours of good quality Y-band spectroscopy during April 
23 (Night 1), and 2.0 hours on April 25 (Night 2 ). 

Data were taken with 180 s exposures and AB dither 
offsets along the slit with ± 1 " and ± 1 . 2 ", respectively. 
In night 2, we also increased the slit width from 077 (as 
used in Night 1) to 079 in anticipation of the slightly 
worse seeing forecast. During these nights we observed 
two masks with a total of eight z ^ 7 — 8 candidate 
galaxies, in addition to lower redshift fillers. 

3.2. Data Reduction 

The data were reduced using a modified version of the 
public MOSEIRE reduction code DRP^. This pipeline 
produces 2D sky-subtracted, rectified, and wavelength- 
calibrated data for each slitlet with a spatial resolution 
of 071799 per pixel and a dispersion of 1.086 A per pixel. 

® https://code.google.com/p/mosfire/ 













Spectroscopic Redshift for a Luminous z = 7.73 Galaxy 


3 




Observed Wavelength |(xm] 


Fig. 2. — Left - Mask layout of the two nights of MOSFIRE 
Y-band observations of our primary target. These two nights pro¬ 
vide two completely independent measurements of this galaxy at 
two different orientations as well as two different positions along 
different slitlets. This also allows us to exclude the possibility of 
contamination in the final stacked spectrum from the two faint 
neighboring galaxies present within 2'' of the primary galaxy along 
the slits. Right - The signal-to-noise ratio around the detected 
emission line in the two independent ID spectra of the two nights, 
averaged over a 4 pixel width (~ 4 A). A line is clearly detected at 
> 4(7 independently in both 2 hr spectra from each night. We also 
checked the unrectified frames to ensure that the positive flux in 
the spectrum indeed originated from the expected position of the 
galaxy along the spectrum. 


Each of our mask contains one slitlct placed on a star 
for monitoring the sky transparency and seeing condi¬ 
tions of each exposure. We use this star to track the mask 
drift across the detector (see, e.g., Kriek et ah 2014), 
which we find to be ±1.5 pixels (±0'.'27) and ±1 pixel 
(±0'.'18) during night 1 and 2, respectively. We sepa¬ 
rately reduce different batches of the data (of 30-45 min 
duration) to limit any S/N reduction caused by this drift, 
before shifting and stacking the data. 

The masks for the two nights have different orienta¬ 
tions (Fig 2). The two independent data sets of the pri¬ 
mary target thus add to the robustness of any detection. 
After creating the 2D spectra for the different masks, 
we applied the appropriate relative shift of the two 2D 
frames before stacking the observations of the two nights 
to our final 2D spectrum. 

Similarly, ID spectra were extracted separately for 
each mask using an optimal extraction based on a pro¬ 
file determined by the slit star. The extracted ID spec¬ 
tra were corrected for Galactic extinction and for telluric 
absorption using nearby AO stars observed in the same 
night at similar airmass. The uncertainty in our opti¬ 
mally extracted ID spectra was determined empirically 
from empty rows in the full, rectified 2D spectra of the 
mask. 

The absolute flux calibration was obtained from the 
slit stars by comparison of the spectra with the 3D-HST 
photometric catalogs (Skelton et al. 2014). An additional 
small correction was applied to account for the extension 
of individual sources in the slit mask by integrating the 
seeing-matched HST images over the slit and comparing 
with the slit loss of stellar sources. 
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f(Lya) = 1.7±0.3 10"''^ erg/s/cm^ 
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Fig. 3.— MOSFIRE spectra of EGS-zs8-l. The full 2D spectrum 
after 2-by-2 binning is shown in the top panel, while the optimally 
extracted ID spectrum is shown on the bottom. The ID spec¬ 
trum was smoothed by a 3 pixel 3 A) moving average filter for 
clarity. The gray shaded area represents the Icr flux uncertainty, 
while the dark red line shows the best-fit model. The line is quite 
extended in the wavelength direction and shows clear asymmetry 
with the expected shape typical for high-redshift Lya lines. The 
spectroscopic redshift measurement is 2;spec = 7.7302 ± 0.0006 in 
excellent agreement with the previously determined photometric 
redshift. Other line characteristics are summarized in Table 1. 


TABLE 1 

Measurements of Galaxy EGS-zs8-1 


Target 

R.A.(J2000) 

14:20:34.89 

Dec (J2000) 

53:00:15.4 

Hieo 

25.03T0.05 

Muv 

-22.06 ± 0.05 

Emission Line 

-S^spec 

7.7302±0.0006 

f(Lya) 

1.7it0.3x 10“^^ erg 

L(LyQ) 

1 . 2 ± 0 . 2 x 10-^3 erg s"! 

EWo(Lya)'‘ 

21±4 A 

Sy, 

15±6 A 

FWHM'^ 

13±3 A 

Vfwhm*’ 

360l7o km s-i 

Physical Parameters'^ 


log Mgai/MQ 9.9T0.2 

logage/yr 8.0±0.5 

log SFR/(Moyr-i) 1.9 ± 0.2 
log SSFR -8.0 ± 0.4 

Auv 1-6 mag 

UV slope {3 —1.7±0.1 


^ Not corrected for IGM absorption. 

^ Derived from truncated Gaussian fit, corrected for 
instrumental broadening, but not for IGM absorp¬ 
tion. 

Based on SED fits (see Sect 5; Oesch et al. 2014). 


4. RESULTS 

Out of the eight z 7 — 8 galaxy candidates, we de¬ 
tected a significant emission line (at > 5cr) for only one 
source (EGS-zs8-l). This line is discussed in detail be¬ 
low. 


4.1. A Lya Emission Line at z = 7.730 

The spectra of our target source EGS-zs8-l (see Ta¬ 
ble 1 for summary of properties) revealed a significant 
emission line at the expected slit position in both masks 
independently (right panels Fig 2). The full 4 hr stacked 
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2D and ID spectra are shown in Figure 3, showing a line 
with a clear asymmetric profile, as expected for a Lya 
line at high redshift ^ 3). Furthermore, it lies at the 
expected wavelength based on our photometric redshift 
estimate 2 ;phot = 7.7 ± 0.3. We therefore interpret this 
line as Lya (other possibilities are discussed in section 
4.3). 

We fit the line using a Markov Chain Monte Carlo 
(MCMC) approach based on the emcee python library 
(Foreman-Mackey et al. 2013). Our model is based on a 
truncated Gaussian profile to account for the IGM ab¬ 
sorption and includes the appropriate instrumental res¬ 
olution. The model also includes the uncertainty on the 
background continuum level. The MGMG output pro¬ 
vides full posterior PDFs and uncertainties for the red¬ 
shift, line flux, significance, and line width. 

The line corresponds to a redshift of ^Lya = 7.7302 ± 
0.0006, with a total luminosity of Li^ya = 1.2 ± 0.2 x 10 ^^ 
erg s“^, and a total detection significance of 6 .la. This 
is somewhat lower, but consistent with a simple estimate 
of 7.2a detection significance from integrating the ID 
extracted pixel flux over the full extent of the line (i.e., 
not accounting for background continuum offsets). 

Note that the redshift of the line is determined from 
our model of a truncated Gaussian profile and is thus 
corrected for instrumental resolution and the asymme¬ 
try arising from the IGM absorption. The peak of the 
observed line (A = 10616 A) thus lies ^ 2.5 A to the red 
of the actual determined redshift. 

Given the brightness of the target galaxy, the de¬ 
tected line corresponds to a rest-frame equivalent width 
EWo = 21 ± 4 A. This is lower than the Lya emitter 
criterion EWq > 25 A set in recent analyses that use the 
Lya fraction among LBGs to constrain the reionization 
process (e.g. Stark et al. 2011; Treu et al. 2013). 

4.2. Line Properties 

Different quantities of the detected line are tabulated 
in Table 1. In particular, we compute the weighted 
skewness parameter, (Kashikawa et al. 2006) find¬ 
ing Sw = 15 ±6 A. This puts the line above the 3 A limit 
found for emission lines at lower redshift (see also section 
4.3). 

The full-width-at-half-maximum (FWHM) of the line 
is quite broad with FWHM = 13 ± 3 A, corresponding 
to a velocity width of Wwhm = 360^70 
galaxy thus lies at the high end of the observed line width 
distribution for 2 ; ^ 5.7 — 6.6 Lya emitters (e.g. Ouchi 
et al. 2010), but is consistent with previous z > 1 Lya 
lines (Ono et al. 2012 ). 

4.3. Caveats 

While the identification of the detected asymmetric 
emission line as Lya is in excellent agreement with the 
expectation from the photometric redshift, we can not 
rule out other potential identifications. As pointed out 
in the previous section, Kashikawa et al. (2006) find that 
weighted asymmetries > 3 A are not seen in lower 
redshift lines, but almost exclusively in Lya of high- 
redshift galaxies. However, at the resolution of our spec¬ 
tra, the observed asymmetry is also consistent with an 
[O II] line doublet in a high electron density environment. 
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Fig. 4. — Top - UV absolute magnitudes of spectroscopically 
confirmed Lyman break galaxies and Lyo; emitters in the cosmic 
reionization epoch, at 2 ; > 6. Our target, EGS-zs8-l (red square), 
represents the highest-redshift source and is the brightest galax¬ 
ies currently confirmed. For reference, the gray dashed line shows 
the evolution of the characteristic magnitude M* of the UV LF 
(Bouwens et al. 2015). The galaxies shown as black squares are as¬ 
sembled from a compilation from Jiang et al. (2013b); Finkelstein 
et al. (2013); Shibuya et al. (2012); Ono et al. (2012), and Vanzella 
et al. ( 2011 ). Bottom - Surface density of the full sample of 2 ; ~ 8 
galaxies in the combined CANDELS and BoRG/HIPPIES fields 
(Bouwens et al. 2015, gray histogram). EGS-zs8-i is the brightest 
and also one of the most massive sources at these redshifts. Note 
that all three 2 : ~ 8 candidates with H ~ 25.0 mag (~ 0.5 mag 
brighter than the rest) are identified in the GANDELS/WIDE sur¬ 
vey area where ancillary V 105 imaging is generally not available. 
Our spectroscopic confirmation is thus especially valuable. 


i.e., with a ratio of [O II]A3726 / [O II]A3729 > 2, and 
with a velocity dispersion of a^ > 100 km s“^. 

If the observed line is an [O II]AA3726, 3729 doublet, 
the redshift of this galaxy would be zqu = 1.85. This 
is very close to the best low redshift SED fit shown in 
Eigure 1. However, that SED requires a strong spectral 
break caused by an old stellar population, for which no 
emission line would be expected. Additionally, the low- 
redshift solution can not explain the extremely red IRAG 
color (used to select this galaxy), and predicts significant 
detections in the ACS/E814W band, as well as in the 
ground-based WIRDS LC-band image (Bielby et al. 2012). 
No such detections are present, however, resulting in a 
likelihood for such an SED of £ < 10“^ (see also Eig 1). 
Thus all the evidence points to this line being Lya at 
^ = 7.73. 


5. DISGUSSION 









Spectroscopic Redshift for a Luminous z = 7.73 Galaxy 


5 


In this Letter we used Keck/MOSFIRE to spectro¬ 
scopically confirm the redshift of one of the brightest 
z 8 galaxies identified by Bouwens et al. (2015) over 
the five CANDELS fields. Interestingly, this source is 
^ 0.5 mag brighter than any source identified in the 
wide-area BoRG and HIPPIES surveys (e.g. Trenti et al. 
2011; Yan et al. 2011; Bradley et al. 2012; Schmidt et al. 
2014). 

As shown in Eig. 4, with ^^gpec = 7.730 and an abso¬ 
lute magnitude Muv = —22.06 ± 0.05 the source EGS- 
zs8-l is currently the most distant and brightest spec¬ 
troscopically confirmed galaxy (apart from a gamma ray 
burst at ^ = 8.2; Tanvir et al. 2009; Salvaterra et al. 
2009). EGS-zs8-l also populates the brightest bin of 
the recent Bouwens et al. (2015) z ^ 8 UV luminosity 
function (LE), which makes it an unusually rare object. 
A spectroscopic confirmation of its high redshift is thus 
particularly valuable for proving the existence of bright 
H = 25.0 mag galaxies at 2 : ^ 8 and for validating the 
bright end LE constraints. 

An SED fit at the spectroscopic redshift of the source 
reveals a relatively high stellar mass log M/Mq = 9.9 ± 
0.2, a star-formation rate of log SER/(MQyr“^) = 1.9 ± 
0.2, and a relatively young, but not extreme age of 
logage/yr = 8.0±0.5 based on an apparent Balmer break 
between the WEC3/IR and the Spitzer photometry (see 
also Table 1). The corresponding formation redshift of 
this galaxy thus lies at Zf = 8.8. Eor details on our SED 
fitting see, e.g., Oesch et al. (2014). 

Interestingly, the source has a UV continuum slope of 
/3 = —1.7±0.1 (measured from the SED fit) and is con¬ 
sistent with considerable dust extinction, E(B—V)= 0.15 
mag. The detection of a significant Lya emission line is 
not inconsistent, however, given the complexities of line 
formation in such young galaxies. 

These observations also allow us to reliably constrain 
the equivalent widths of the [O IIl]±H/d emission lines 
in this galaxy based on its IRAC colors. At the spec¬ 
troscopic redshift of the source, these lines are shifted in 
the 4.5 yum channel resulting in a color of [3.6] —[4.5] = 
0.53 ±0.09 mag. This is consistent with a combined rest- 
frame equivalent width of EWo([0 IIl]±H/3)= 7201^50 

A. 

Such a high equivalent width is in some tension 
with the inferred stellar population age of ^ 100 Myr. 
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However, it is completely consistent with the average 
EWo([0 IIl]±H;d) found for z ^ 7 — 8 galaxies in previ¬ 
ous work (e.g. Labbe et al. 2013; Smit et al. 2014, 2015), 
where such strong lines were found to be ubiquitous (see 
also Laporte et al. 2014). 

The fact that these strong lines are seen in a significant 
fraction of the z 7 — 8 galaxies is at odds with the 
interpretation of extremely young galaxy ages of < 10 
Myr (e.g. Einkelstein et al. 2013). Nevertheless, very 
stochastic star-formation may explain some of the line 
strength. Instead, it is possible that an evolution in the 
ionization properties of early galaxy populations may be 
causing stronger emission lines with more extreme line 
ratios [O IIl]/H/3 as is observed at 2 ; ^ 3 — 4 (see e.g. 
Holden et al. 2014). We will discuss models that will 
provide greater insights into these strong emission line 
sources in a future paper. 

Our confirmation of a source with extremely strong 
rest-frame optical emission lines at 2:spec = 7.730 together 
with two very similar sources at ^^gpec = 7.213 and 7.508 
(Ono et al. 2012; Einkelstein et al. 2013) provides fur¬ 
ther support for the likelihood of ubiquitous strong rest- 
frame optical lines as evidenced in the IRAG photometry 
of ^ ^ 7 — 8 galaxies. This has important consequences 
for future observations with JWST, which promises ex¬ 
tremely efficient spectroscopic follow-up of such strong 
line emitters with NIRspec out to the highest redshifts 
of currently known galaxies. 
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